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REMARKS 

The specification has been amended to correct typographical errors. Applicants submit 
that no new matter was added by the amendments. 

Claim 27 was amended to remove reference to the figure. Claims 27 and 32-35 remain 
pending. Applicants thank the Examiner for her careful review of the instant application and for 
withdrawing the rejections over the prior art. The pending claims stand rejected for allegedly 
being unpatentable under 35 U.S.C. §§101 and 112, first paragraph. For the reasons set forth 
below, Applicants respectfully traverse. 

Rejection under 35 U.S.C. §101 

The PTO maintains the rejection of all claims under 35 U.S.C. § 101 as lacking a specific, 
substantial, and credible utility. The PTO does not accept the assertion that the gene 
amplification data of Example 16 establish that the PRO 1800 polypeptide is overexpressed in 
certain cancers and is therefore useful as a cancer diagnostic tool. 

Applicants respectfully traverse. 
Utility - Evidentiary Standard 

An Applicant's assertion of utility creates a presumption of utility that will be sufficient to 
satisfy the utility requirement of 35 U.S.C. § 101, "unless there is a reason for one skilled in the 
art to question the objective truth of the statement of utility or its scope." In re Langer, 503 F.2d 
1380, 1391, 183 USPQ 288, 297 (CCPA 1974). See, also In reJolles, 628 F.2d 1322, 206 USPQ 
885 (CCPA 1980); In re Irons, 340 F.2d 974, 144 USPQ 351 (1965); In re Sichert, 566 F.2d 
1154, 1159, 196 USPQ 209, 212-13 (CCPA 1977). 

Compliance with 35 U.S.C. § 101 is a question of fact. Raytheon v. Roper, 724 F.2d 951, 
956, 220 USPQ 592, 596 (Fed. Cir. 1983) cert, denied, 469 US 835 (1984). According to the 
PTO Utility Examination Guidelines (2001), irrefutable proof of a claimed utility is not required. 
Rather, a specific, substantial, and credible utility requires only a "reasonable" confirmation of a 
real world context of use. The evidentiary standard to be used throughout ex parte examination 
in setting forth a rejection is a preponderance of the totality of the evidence under consideration. 
In re Oetiker, 977 F.2d 1443, 1445, 24 USPQ2d 1443, 1444 (Fed. Cir. 1992). The M.P.E.P. 
§ 2107 clearly state that "[w]here the asserted utility is not specific or substantial, a prima facie 



-7- 



Appl. No. 
Filed 



09/866,034 
May 25, 2001 



showing must establish that it is more likely than not that a person of ordinary skill in the art 
would not consider that any utility asserted by the applicant would be specific and substantial." 
(emphasis added). The M.P.E.P. also states that the prima facie showing must contain the 
"[sjupport for the factual findings relied upon" in reaching the conclusion that the utility is not 
substantial and specific. Thus, the initial burden is on the PTO to establish that it is more likely 
than not that one of skill in the art would not consider the asserted utility substantial, and that 
conclusion must be supported. Only after the PTO has made a properly supported prima facie 
showing of lack of utility does the burden of rebuttal shift to the applicant. The issue will then be 
decided on the totality of evidence, with the applicant needing to provide a "reasonable" 
confirmation of a real world context of use by a preponderance of the totality of the evidence. 
The Art supports the Applicants ' Assertion that Gene Amplification leads to Increased Gene 
Expression 

Applicants first address the PTO's conclusion that the ACt data does not establish that the 
gene encoding PRO 1800 is overexpressed. Rather, the PTO asserts that the data show, at best, 
that the nucleic acid is present at levels only twice that of the control cells. The PTO asserts that 
this is not conclusive or even indicative that the gene is overexpressed. Instead, the PTO 
suggests that the more likely explanation is that the two-fold increase in the gene is indicative of 
aneuploidy, not overexpression. The PTO asserts that if the result is due to aneuploidy, then a 
person of ordinary skill in the art would expect the level of PRO 1800 polypeptide, if present at 
all, to not be affected at all, or only in a minor fashion. 

Applicants submit that those of skill in the art would recognize that amplification of a 
gene due to aneuploidy will more likely than not lead to an increase in expression of that gene as 
measured by the level of mRNA. This assertion is supported by numerous references. Orntoft et 
al. {Molecular and Cellular Proteomics, 1:37-45 (2002)) (submitted herewith as Exhibit 1) 
studied transcript levels of 5600 genes in malignant bladder cancers which were linked to a 
gain/loss of chromosomal material using an array-based method. Orntoft et al showed that there 
was a gene dosage effect and teach that "in general (18 of 23 cases) chromosomal areas with 
more than 2-fold gain of DNA showed a corresponding increase in mRNA transcripts" (Orntoft 
at 37, column 1, abstract). In addition, Hyman et al {Cancer Research, 62:6240-6245 (2002)) 
(submitted herewith as Exhibit 2) used CGH analysis and cDNA microarrays to compare DNA 
copy numbers and mRNA expression of over 12,000 genes in breast cancer tumors and cell lines. 
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They showed that there is "evidence of a prominent global influence of copy number changes on 
gene expression levels" (Hyman at 6244, column 1, last paragraph). 

Additional supportive teachings are also provided by Pollack et al. (PNAS, 99:12963- 
12968 (2002)) (submitted herewith as Exhibit 3) who studied a series of primary human breast 
tumors and found that "[b]y analyzing mRNA levels in parallel, we have also discovered that 
changes in DNA copy number have a large, pervasive, direct effect on global gene expression 
patterns in both breast cancer cell lines and tumors." (Pollack at 12967 at column 1, emphasis 
added). Their study found that "62% of highly amplified genes show moderately or highly 
elevated expression, that DNA copy number influences gene expression across a wide range of 
DNA copy number alterations (deletion, low-, mid- and high-level amplification), that on 
average, a 2-fold change in DNA copy number is associated with a corresponding 1.5 -fold 
change in mRNA levels." (Pollack at 12963, column 1, abstract). This report is particularly 
persuasive because the high-resolution comparative genomic hybridization analysis used to 
assess DNA copy number was particularly sensitive. 

Together, these articles collectively teach that in general, gene amplification increases 
mRNA expression. They directly support the Applicants' assertion that the amplification of the 
PRO 1800 gene reported in Example 16, even where the amplification is only two-fold, leads to 
an increase in expression of the PRO 1800 gene. 

The Pennica reference (PNAS USA, 95:14717-14722 (1998)) cited by the PTO does not 
support the PTO's argument that because gene amplification is not necessarily correlated to 
mRNA level, there is no correlation between nucleic acid levels and protein expression. The 
PTO focuses on the statement from Pennica that the WISP-2 gene DNA was amplified in colon 
tumors, but mRNA expression was reduced in the majority of the tumors. (Pennica at 14722, 
column 1). This inverse correlation is in contrast to the WISP-1 gene, which was amplified and 
had higher mRNA levels. The authors of Pennica offer an explanation for what they obviously 
viewed as an anomalous result: "Because the center of the 20ql3 amplicon [of which WISP-2 is a 
part] has not yet been identified, it is possible that the apparent amplification observed for WISP- 
2 may be caused by another gene in this amplicon." (Pennica at 14722, column 1). Thus, the 
example of a lack of positive correlation between gene amplification and RNA levels relied on 
by the PTO may not actually be real. The fact that the authors attempt to explain this anomaly 
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only supports Applicants' argument that the accepted understanding in the art is that there is a 
direct correlation between gene amplification and the level gene expression. 

Even if the W1SP-2 result is a real example of gene amplification not leading to increased 
gene expression, such a result only leads to the conclusion that increased gene copy number does 
not necessarily result in increased gene expression. The standard for establishing utility, 
however, is not absolute certainty, but rather that it is more likely than not that one of skill in the 
art would believe the asserted utility. The fact that in one case there seemed to be no correlation 
between gene amplification and the level of mRNA expression does not mean that the PTO has 
established that it is more likely than not that, in general, such correlation does not exist. The 
PTO has not shown whether the lack of correlation observed for one member of the family of 
WISP polypeptides is typical, or is merely a discrepancy, an exception to the rule of correlation. 
Indeed, the authors' attempt to explain the lack of correlation suggests that this result is the 
exception, and not the rule. 

Given the teaching of the references submitted herewith, Applicants submit that they have 
established that one of skill in the art would believe that it is more likely than not that the 
reported amplification of the PRO 1800 gene will lead to increased gene expression. As 
Applicants demonstrate below, it is well-established in the art that there is a direct correlation 
between gene expression levels and the level of the encoded protein. When combined, these 
teachings establish a "reasonable" confirmation that the claimed polypeptides are overexpressed 
in certain cancers, and are therefore useful as cancer diagnostic tools. 

A pplicants have established that the Accepted Understanding in the Art is that there is a Direct 
Correlation between mRNA Levels and the Level of Expression of the Encoded Protein 

The Applicants have previously offered the Declaration of J. Christopher Grimaldi, an 
expert in the field, to support their assertion that it is well-established in the art that there is a 
direct correlation between increased gene expression and increased expression of the encoded 
protein. The PTO has rejected the assertion in paragraph 5 of the Grimaldi declaration that 
increased mRNA expression is expected to be associated with increased protein production, 
stating that it appears to be the declarant's opinion and is not supported by fact or evidence. The 
PTO also asserts that the specification and the record in general do not distinguish between total 
nucleic acid (which includes chromosomal DNA) and mRNA. The PTO states that one cannot 
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determine from the specification if the amplification reported in Example 16 is due to an increase 
in chromosomal copy number or increased transcription rates. 

First, Applicants submit that it is clear from the Specification that the amplification 
reported in Example 16 is a measure of chromosomal DNA, and is due to an increase in copy 
number. Applicants apologize for any confusion in the record caused by any failure to 
distinguish between chromosomal DNA and mRNA. 

Next, Applicants respectfully submit that the previously entered Declaration of J. 

Christopher Grimaldi supports the Applicants position that it is well-established that there is a 

direct correlation between increased gene expression and expression of the encoded protein. In 

paragraph 5 of his declaration, Grimaldi states that: 

Those who work in this field are well aware that in the vast majority of cases, 
when a gene is over-expressed,... the gene product or polypeptide will also be 
over-expressed. . . . Techniques used to detect mRNA, such as Northern Blotting, 
Differential Display, in situ hybridization, quantitative PCR, Taqman, and more 
recently Microarray technology all rely on the dogma that a change in mRNA will 
represent a similar change in protein. If this dogma did not hold true then these 
techniques would have little value and would not be so widely used. 

Applicant reminds the Examiner that "Office personnel must accept an opinion from a 
qualified expert that is based upon relevant facts whose accuracy is not being questioned; it is 
improper to disregard the opinion solely because of a disagreement over the significance or 
meaning of the facts offered." PTO Utility Examination Guidelines (2001). Grimaldi's 
conclusion that there is a general correlation between mRNA and protein levels is supported by 
the fact that all of the techniques used to measure mRNA levels would have little utility if such a 
correlation did not exist. 

Applicants also submit herewith (as Exhibit 4) a copy of the Declaration of Dr. Paul 

Polakis, an expert in the field of cancer biology and principal investigator of the Tumor Antigen 

Project of Genentech, Inc., the assignee of the present application, originally submitted in a 

related and co-owned patent application Serial No. 10/032,996. As Dr. Polakis explains, the 

primary focus of the microarray project was to identify tumor cell markers useful as targets for 

both the diagnosis and treatment of cancer in humans. The scientists working on the project 

extensively rely on results of microarray experiments in their effort to identify such markers. As 

Dr. Polakis explains, using microarray analysis, Genentech scientists have identified 

approximately 200 gene transcripts (mRNAs) that are present in human tumor cells at 
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significantly higher levels than in corresponding normal human cells. To date, they have 
generated antibodies that bind to about 30 of the tumor antigen proteins expressed from these 
differentially expressed gene transcripts and have used these antibodies to quantitatively 
determine the level of production of these tumor antigen proteins in both human cancer cells and 
corresponding normal cells. Having compared the levels of mRNA and protein in both the tumor 
and normal cells analyzed, they found a very good correlation between mRNA and corresponding 
protein levels. Specifically, in approximately 80% of their observations they have found that 
increases in the level of a particular mRNA correlates with changes in the level of protein 
expressed from that mRNA. 

While the proper legal standard is to show that the existence of correlation between 
mRNA and polypeptide levels is more likely than not, the showing of approximately 80% 
correlation for the molecules tested in the Polakis Declaration greatly exceeds this legal standard. 

As stated in paragraph 6 of his declaration: 

Based on my own experience accumulated in more than 20 years of research, 
including the data discussed in paragraphs 4 and 5 above [showing a positive 
correlation between mRNA levels and encoded protein levels in the vast majority 
of cases] and my knowledge of the relevant scientific literature, it is my 
considered scientific opinion that for human genes, an increased level of mRNA 
in a tumor cell relative to a normal cell typically correlates to a similar increase in 
abundance of the encoded protein in the tumor cell relative to the normal cell. In 
fact, it remains a central dogma in molecular biology that increased mRNA levels 
are predictive of corresponding increased levels of the encoded protein. 
(Emphasis added). 

Applicants submit that Dr. Polakis' opinion that it is a "central dogma in molecular biology that 
increased mRNA levels are predictive of corresponding increased levels of the encoded protein" 
is well-supported by the facts outlined above, and therefore should not be disregarded by the 
PTO. 

Finally, Applicants note that the above statements from the declarations of Grimaldi and 
Polakis are supported by the teaching of Orntoft et al In addition to determining that there is a 
correlation between gene amplification and increased gene expression, Orntoft et al also 
examined the correlation between gene expression and protein expression. They stated that 
"[w]ith few exceptions we found a good correlation (p<0.005) between transcript alteration and 
protein levels." (Orntoft at 37, column 1, abstract). 
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Taken together, this reference and the Grimaldi and Polakis declarations establish that the 
accepted understanding in the art is that an increase in the level of mRNA leads to an increase in 
the level of the encoded protein. Therefore, Applicants submit that they have established that it 
is more likely than not that one of skill in the art would believe that the PRO 1800 polypeptide is 
overexpressed in certain cancers. 

The combined teachings in the art, exemplified by Orntoft et al, Hyman et al. and 
Pollack et al and the Grimaldi and Polakis declarations, overwhelmingly teach that gene 
amplification influences gene expression and that gene expression influences protein levels. 
Thus, one of skill in the art would reasonably expect, in this instance, based on the gene 
amplification data for the PRO 1800 gene, that the PRO 180 protein is concomitantly 
overexpressed in certain cancers. Thus, Applicants submit that they have provided "reasonable" 
confirmation that the PRO 1800 proteins have utility in the diagnosis of cancer. 
The Claimed Polypeptides would have Diagnostic Utility even if there is no Direct Correlation 
between Gene Amplification and Protein Expression 

Applicants believe they have supplied sufficient evidence to show that there is a 
significant correlation between gene amplification and protein expression. Even assuming 
arguendo that, there is no direct correlation between gene amplification and protein expression 
for PRO 1800, which Applicants submit is not true, a polypeptide encoded by a gene that is 
amplified in cancer would still have a credible, specific and substantial utility. 

As explained in paragraph 6 of the previously submitted Ashkenazi Declaration: 

Even when amplification of a cancer marker gene does not result in significant 
over-expression of the corresponding gene product, this very absence of gene 
product over-expression still provides significant information for cancer 
diagnosis and treatment. Thus, if over-expression of the gene product does not 
parallel gene amplification in certain tumor types but does so in others, then 
parallel monitoring of gene amplification and gene product over-expression 
enables more accurate tumor classification and hence better determination of 
suitable therapy. In addition, absence of over-expression is crucial information 
for the practicing clinician. If a gene is amplified but the corresponding gene 
product is not over-expressed, the clinician accordingly will decide not to treat a 
patient with agents that target that gene product. 

This statement is echoed by Grimaldi in his declaration at paragraph 6, and is further 
supported by the teachings in the article by Hanna and Mornin (submitted herewith as Exhibit 5). 
The article teaches that the HER-2/neu gene has been shown to be amplified and/or 
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overexpressed in 10%-30% of invasive breast cancers and in 40-60% of intraductal breast 
carcinoma. Further, the article teaches that diagnosis of breast cancer includes testing both the 
amplification of the HER-2/neu gene (by FISH) as well as the overexpression of the HER-2/neu 
gene product (by IHC). Even when the protein is not overexpressed, the assay relying on both 
tests leads to a more accurate classification of the cancer and a more effective treatment of it. 

The Applicants have established that it is the general, accepted understanding in the art 
that gene amplification leads to increased gene and protein expression. Because this is the 
general rule, even when this is not the case, a polypeptide encoded by a gene that is 
overexpressed in cancer would still have utility. Thus, Applicants have demonstrated another 
"reasonable" confirmation of the asserted utility for the claimed polypeptides. 

Conclusion 

Applicants have provided several references and expert opinions supporting the utility of 
the present invention. Applicants have established, supported by three references, that the 
general rule accepted in the art is that gene amplification leads to an increase in gene expression. 
Thus, Applicants have provided sufficient evidence that one of skill in the art is more likely than 
not to accept that the evidence of amplification of the PRO 1800 gene is also evidence of 
overexpression of the gene. Similarly, Applicants have offered the opinions of two experts in the 
field, each supported by relevant facts, which state that increased gene expression is understood 
by those of skill in the art to lead to increased protein expression for the encoded protein. These 
opinions are supported by an additional reference. Thus, Applicants have provided sufficient 
evidence that one of skill in the art is more likely than not to accept that the evidence of 
amplification of the PRO 1800 gene and the accompanying increase in gene expression is 
evidence that the PRO 1800 protein is overexpressed in certain cancers. Applicants submit that it 
is beyond dispute that a protein which is overexpressed in cancer has utility as a cancer 
diagnostic tool, as well as a target for therapy. 

Finally, even if one does not accept that gene amplification leads to increased protein 
expression, Applicants have submitted the opinion of two experts in the field, along with 
supporting art, that those of skill in the art would still find it useful to monitor the level of the 
gene amplified in cancer, as well as the level of the encoded protein. This is done to further 
classify the type of cancer, as well as select the most effective treatment. Thus, even if the 
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amplification of the PRO 1800 gene in certain cancers does not lead to overexpression of the 
PRO 1800 protein, the claimed polypeptides would still have utility as a cancer diagnostic tool. 

Rejection under 35 U.S.C. SI 12, first paragraph 

Claims 27, 28, and 32-35 were also rejected under 35 U.S.C. §1 12, first paragraph as not 
being enabled because they are not supported by an asserted utility, and therefore one would not 
know how to use the claimed invention. For the reasons outlined above in response to the 
rejection under 35 U.S.C. §101, Applicants believe they have established a utility of the claimed 
invention and therefore respectfully request withdrawal of the rejection under 35 U.S.C. §112, 
first paragraph. 



In view of the above, Applicants respectfully maintain that claims are patentable and 
request that they be passed to issue. Applicants invite the Examiner to call the undersigned if any 
remaining issues may be resolved by telephone. 

Please charge any necessary fees, including any fees for any extensions of time, to 
Deposit Account No. 11-1410. 



CONCLUSION 



Respectfully submitted, 



KNOBBE, MARTENS, OLSON & BEAR, LLP 
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Genome-wide Study of Gene Copy Numbers, 
Transcripts, and Protein Levels in Pairs of 
Non-invasive and Invasive Human Transitional 
Cell Carcinomas* 

Torben F. 0rntoft|:§ f Thomas ThykJaetU, Frederic M. Waldman||, Hans Wolf**, 
and Julio E. Cells** 



Gain and loss of chromosomal material is characteristic 
of bladder cancer, as well as malignant transformation in 
general. The consequences of these changes at both the 
transcription and translation levels Is at present unknown 
partly because of technical limitations. Here we have at- 
tempted to address this question in pairs of non-invasive 
and Invasive human bladder tumors using a combination 
of technology that included comparative genomic hybrid- 
ization, high density oligonucleotide array-based monitor- 
ing of transcript levels (5600 genes), and high resolution 



phenomenon at both the transcription and translation levels. 
High throughput array studies of the breast cancer cell line 
BT474 has suggested that there is a correlation between 
DNA copy numbers and gene expression in highly amplified 
areas (2), and studies of individual genes in solid tumors 
have revealed a good correlation between gene dose and 
mRNA or protein levels in the case of c-erb-B2, cyclft) d1 t 
emsl, and N-myc (3-5). However, a high cyclin D1 protein 
expression has been observed without simultaneous am- 



two-dimensional gel electrophoresis/the results showed z^P^ 03 ** 0 " W). and a low level of c-myc copy number in 



that there is a gene dosage effect ihat in some cases 
superimposes on other regulatory mechanisms. This ef- 
fect depended (p < 0.01?) on the magnitude of the com- 
parative genomic hybridization change. In general (18 of 
23 cases), chromosomal areas with more than 2-foid gain 
of DNA showed a corresponding increase in mRNA tran- 
scripts. Areas with loss of DNA, on the other hand, 
showed either reduced or unaltered transcript levels) Be- 
cause most proteins resolved by two-dimensional gets 
are unknown it was only possible to compare mRNA and 
protein alterations ta relatively few cases of well focused 
abundant proteins. VVith few exceptions we found a good 
correlation (p < 0*005) between transcript alterations and 
protein levels. The implications, as well as limitations, 
of the approach are discussed. Molecular & Cellular 
Proteomlcs 1:37-45, 2002. 



, Aneuploidy is a common feature of most human cancers 
(1), but little is known about the genome-wide effect of this 
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crease was observed without concomitant c-myc protein 
overexpression (6). 

In human bladder tumors, karyotyping, fluorescent In situ 
hybridization, and comparative genomic hybridization (CGH) 1 
have revealed chromosomal aberrations that seem to be 
characteristic of certain stages of disease progression. In the 
case of non-invasive pTa transitional ceil carcinomas (TCCs), 
this includes loss of chromosome 9 or parts of it, as well as 
loss of Y in males. In minimally invasive pT1 TCCs, the fol- 
lowing alterations have been reported: 2q-, 11p», 1q+, 
11q13+ t 17q+, an^J 20q+ (7-12). It has been suggested that 
these regions harbor tumor suppressor genes and onco- 
genes; however, the large chromosomal areas Involved often 
contain many genes, making meaningful predictions of the 
functional consequences of losses and gains very difficult. 

In this investigation we have combined gehome-wide tech- 
nology for detecting genomic gains and losses (CGH) with 
gene expression profiling techniques (microarrays and pro- 
teomlcs) to determine the effect of gene copy number on 
transcript and protein levels in pairs of non-Invasive and In- 
vasive human bladder TCCs. 

EXPERIMENTAL PROCEDURES 

Materia/— Bladder tumor biopsies were sampled after Informed 
consent was obtained and after removal of tissue for routine pathol- 
ogy examination. By light microscopy tumors 335 and 532 were 
staged by an experienced pathologist as pTa (superficial papillary), 

1 The abbreviations used are: CGH, comparative genomic hybrid- 
ization; TCC, transitional cell carcinoma; LOH, loss of heterozygosity; 
PA-FABP, psoriasis-associated fatty acid-binding protein; 2D, 
two-dimensional. 
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L °° Py number 0X16 mRNA expression level. Shown from left to right are chromosome (Cnr.), GGH profiles, gene location and 
expression level of specific genes, and overall expression level along the chromosome- A, expression of mRNA In Invasive tumor 733 as 
compared with the non-invasive counterpart tumor 335. fl, expression of mRNA in invasive tumor 827 compared with the non-invasive 
c^merpart tumor 532. The average fluorescent signal ratio between tumor DMA and normal DNA is shown along the length of the chromosome 
Vm. The bold carve in the ratio profile represents a mean of four chromosomes and is surrounded by thin curves indicating one standard 
n c V » Hon* The central vertical fae broken) indicates a ratio value of 1 (no change), and the vertical lines next to It (dotted) indicate a ratio of 
0.5 (/eft) and 2.0 (right). In chromosomes where the non-Invasive tumor 335 used for comparison showed alterations In DNA content, the ratio 
profile of that chromosome is shown to the right of the invasive tumor profile. The colored bars represents one gene each, Identified by the 
running numbers above the bars (the name of the gene can be seen at www.MDLDK/sdata.html). The bars Indicate the purported location of 
the gene, and the colors Indicate the expression level of the gene in the invasive tumor compared with the non-Invasive counterpart; >2-fold 
tocrease (WacA), >2-fo!d decrease (Wue) f no significant change {orange). The oar to the far right, entitled Expression shows the resulting change 
m expression along the chromosome; the colors Indicate that at least half of the genes were up-regulated (Mac*), at least half of the genes 
down-f eguteted (blue), or more than half of the genes are unchanged {orange). If a gene was absent In one of the samples and present In 
another, it was regarded as more than a 2-fold change. A 2-fold level was chosen as this corresponded to one standard deviation In a double 
determination of -1800 genes. Centromeres and heterochromatic regions were excluded from data analysis. 



Grade I and ll, respectively, tumors 733 and 827 were staged as pT1 
(invasive Into submucosal 733 was staged as solid, and 827 was 
staged as papillary, both grade III. 

mRNA Preparation ~~l}ssue biopsies, obtained fresh from surgery, 
were embedded immediately In a sodium-guanldinlum thfocyanate 
solution and stored at -80 °C. Total RNA was Isolated using the 
RNAzol B RNA isolation method (WAK-Chemie Medical GMBH). 
poly(A) + RNA was Isolated by an oligo(dT) selection step (Ollgotex 
mRNA kit; Qiagen). 

cRNA Preparation-! ^g of mRNA was used as starting material. 
The first and second strand cDNA synthesis was performed using the 
Superscript® choice system (Invitrogen) according to the manufac- 
turer's Instructions but using an oligo(d7) primer containing a T7 RNA 
polymerase binding site. Labeled cRNA was prepared using the ME- 
GAscrip® in vitro transcription kit (Ambion). Biotin-labeled CTP and 



UTP (Enzo) was used, together with unlabeled NTPs in the reaction. 
Following the in vitro transcription reaction, the unincorporated nu- 
cleotides were removed using RNeasy columns (Qiagen). 

Array Hybridization and Scanning— Array hybridization and scan- 
ning was modified from a previous method (13). 10 j*g of cRNA was 
fragmented at 94 °C for 35 min In buffer containing 40 mM Tris 
acetate, pH 8.1, 100 mM KOAc, 30 mM MgOAc. Prior to hybridization, 
the fragmented cRNA in a 6x SSPE-T hybridization buffer (1 m NaCI, 
10 mM Tris, pH 7.6, 0.005% Triton), was heated to S5 °C for 5 min, 
subsequently cooled to 40 "0, and loaded onto the Affymetrix probe 
array cartridge. The probe array was then incubated for 16 h at 40 °C 
at constant rotation (60 rpm). The probe array was exposed to 10 
washes In 6X SSPE-T at 25 °C followed by 4 washes in 0.5X SSPE-T 
at 50 *C. Tne biotinylated cRNA was, stained, with a streptavldin- 
phycoerythrin conjugate, 10 jug/ml (Molecular Probes) in 6x SSPE-T 
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Fig. 1 — continued 



for 30 min at 25 °C followed by 1 0 washes in 6x SSPE-T at 25 ft C. The 
probe arrays were scanned at 560 nm using a confocal laser scanning 
microscope (made for Affymetrtx by Hewlett-Packard). The readings 
from the quantitative scanning were analyzed by Affymetrix gene 
expression analysis software. 

MiavsateWte Anafysfc— MicrosateRite Analysis was performed as 
described previously (14). Microsateliites were selected by use of 
wwwJKbLnlm.nih.gdv/genemap98, and primer sequences were ob- 
tained from the genome data base at www.gcfo.org. DNA was extracted 
from tumor and blood and amplified by PCR In a volume of 20 /d for 35 
cycles. The ampilcons were denatured and electrophoresed for 3 h In an 
ABi Prism 377. Data were collected in the Gene Scan program for 
fragment analysis. Loss of heterozygosity was defined as less than 33% 
of one allele detected in tumor amplicons compared with blood, 

Proteomic Ana!ysis-TCCs were minced Into small pieces and 
homogenized in a small glass homogenizer In 0*5 ml of lysis solution- 
Samples were stored at -20 °C until use. The procedure for 2D gel 
electrophoresis has been described in detail elsewhere (15, 16). Gels 
were stained with silver nitrate and/or Coorhassie Brilliant Blue. Pro- 
teins were Identified by a combination of procedures that Included 
microsequencing, mass spectrometry, two-dimensional gel Western 
immunoblotting, and comparison with the master two-dimensional gel 
Image of human keratinocyte proteins; see biobase.dk/cgi-bin/cefis. 

CGH— Hybridization of differentially labeled tumor and normal DNA 
to normal metaphase chromosomes was performed as described 
previously (10), Fluorescein-labeted tumor DNA (200 ng), Texas Red- 



labeled reference DNA (200 ng), and human Cot-1 DNA (20 ^g) were 
denatured at 37 °C for 5 min and applied to denatured normal met- 
aphase slides. Hybridization was at 37 °C for 2 days. After washing, 
the slides were counterstained with 0.1 5 ftg/ml 4 f 0-dtamldino-2-phe- 
nytindote in an anti-fade solution. A second hybridization was per- 
formed for alt tumor samples using fluoresceJn-tabeled reference DNA 
and Texas Red-labeled tumor DNA (inverse labeling) to confirm the 
aberrations detected during the initial hybridization. Each CGH ex- 
periment also included a normal control hybridization using fluores- 
cein- and Texas Red-labeled normal DNA. Digital image analysis was 
used to identify chromosomal regions with abnormal fluorescence 
ratios, Indicating regions of DNA gains and losses. The average 
greercred fluorescence intensity ratio profiles were calculated using 
four images of each chromosome (eight chromosomes total) with 
normalization of the green:red fluorescence Intensity ratio for the 
entire metaphase and background correction Chromosome identifi- 
cation was performed based on 4,6-dlamldlno-2-phenylindole band- 
ing patterns. Only images showing uniform high intensity fluores- 
cence with minimal background staining were analyzed. All 
centromeres, p arms of acrocentric chromosomes, and heterochro- 
matic regions were excluded from the analysis, 

RESULTS 

Comparative Genomic Hybridization— The CGH analysis 
identified a number, of chromosomal gains and losses in the 
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Table I 

Correlation between alterations detected by CGH and by expression monitoring 

Top, CGH used as independent variable (if CGH alteration - what expression ratio was found); bottom, altered expression used as 
independent variable (if expression alteration - what CGH deviation was found). 
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two Invasive tumors (stage pT1 , TCCs 733 and 827), whereas 
the two non-invasive papillomas {stage pTa, TCCs 335 and 
532) showed only 9p-, 9q22-q33-, and X-, and 7+, 9q-, 
and Y-, respectively. Both invasive tumors showed changes 
(1q22-24+, 2q14.1-qter- t 3q12-q13.3- 6q12-q22~, 
9q344 , 11q12-q13+, 17+, and 20q11.2-q12+) that are typ- 
ical for their disease stage, as well as additional alterations, 
some of which are shown in Fig. 1. Areas with gains and 
losses deviated from the normal copy number to some extent, 
and the average numerical deviation from normal was 0.4-fold 
in the case of TCC 733 and 0.3-fold for TCC 827. The largest 
changes, amounting to at least a doubling of chromosomal 
content were observed at 1q23 In TCC 733 (Fig. 1A) and 
20q12 in TCC 827 (Fig. 16). 

mRNA Expression in Relation to DNA Copy Number— The 
mRNA levels from the two invasive tumors (TCCs 827 and 
733) were compared with the two non-invasive counterparts 
(TCCs 532 and 335). This was done in two separate experi- 
ments in which we compared TCCs 733 to 335 and 827 to 
532, respectively, using two different scaling settings for the 
arrays to rule out scaling as a confounding parameter. Ap- 
proximately 1,800 genes that yielded a signal on the arrays 
were searched in the Unigene and Genemap data bases for 
cfcronTosomaJ location, and those with a known location 
(1096) were plotted as bars covering their purported locus. In 
that way It was possible to construct a graphic presentation of 
DNA copy number and relative mRNA levels along the indi- 
vidual chromosomes (Fig. 1). 

For each mRNA a ratio was calculated between the level in 
the invasive versus the non-invasive counterpart Bars, which 
represent chromosomal location of a gene, were color-coded 
according to the expression ratio, and only differences larger 



than 2-fold were regarded as informative (Fig. 1). The density 
of genes along the chromosomes varied, and areas contain- 
ing only one gene were excluded from the calculations. The 
resolution of the CGH method is very tow, and some of the 
outlier data may be because of the feet that the boundaries of 
the chromosomal aberrations are not known at high resolution. 

Two sets of calculations were made from the data. For the 
first set we used CGH alterations as the independent variable 
and estimated the frequency of expression alterations in these 
chromosomal areas. In general, areas with a strong gain of 
chromosomal material contained a cluster of genes having 
increased mRNA expression. For example, both chromo- 
somes 1q21-q25, 2p and 9q, showed a relative gain of more 
than 100% in DNA copy number that was accompanied by 
increased mRNA expression levels in the two tumor pairs (Fig. 
1). In most cases, chromosomal gains detected by CGH were 
accompanied by an increased level of transcripts in both 
TCCs 733 (77%) and 827 (80%) (Table I, fop). Chromosomal 
losses, on the other hand, were not accompanied by de- 
creased expression in several cases, and were often regis- 
tered as having unaltered RNA levels (Table I, top). The inabil- 
ity to detect RNA expression changes in these cases was not 
because of fewer genes mapping to the lost regions (data not 
shown). 

In the second set of calculations we selected expression 
alterations above 2-fold as the Independent variable and es- 
timated the frequency of CGH alterations In these areas. As 
above, we found that increased transcript expression corre- 
lated with gain of chromosomal material (TCC 733, 69% and 
TCC 827, 59%), whereas reduced expression was often de- 
tected in areas with unaltered CGH ratios (Table I, bottom). 
Furthermore, as a control we looked at areas with no alter- 
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Fig. 2. Correlation between maximum CGH aberration and the abifity to detect expression change by oligonucleotide array 
monitoring. The aberration is shown as a numerical -fold change in ratio between invasive tumors 827 (A) and 733 (♦) and their non-invasive 
counterparts 532 and 335. The expression change was taken from the Expression line to the right In Rg. 1, which depicts the resulting 
expression change for a given chromosomal region. At least half of the mRNAs from a given region have to be either up- or down-regulated 
to be scored as an expression change. All chromosomal arms in which the CGH ratio plus or minus one standard deviation was outside the 
ratio value of one were included. 



ation in expression. No alteration was detected by CGH in 
most of these areas (TCC 733, 60% and TCC 827, 81 %; see 
Table I, bottom). Because the ability to observe reduced or 
Increased mRNA expression clustering to a certain chromo- 
somal area clearly reflected the extent of copy number 
changes, we plotted the maximum CGH aberrations In the 
regions showing CGH changes against the ability to detect a 
change in mRNA expression as monitored by the oligonucleo- 
tide arrays (Rg, 2)(&)r both tumors TCC 733 (p < 0.01 5) and 
TCC 827 (p < 0.00003) a highly significant correlation was 
observed between the level of CGH ratio change (reflecting 
. the DNA copy number) and alterations detected by the array 
based technology (Fig. 2| Similar data were obtained when 
areas with altered expression were used as Independent vari- 
ables. These areas correlated best with CGH when the CGH 
ratio deviated 1.6- to 2-0-fold (Table I, bottom) but mostly did 
not at lower CGH deviations. These data probably reflect that 
loss of an allele may only lead to a 50% reduction In expres- 
sion level, which is at the cut-off point for detection of expres- 
sion alterations. Gain of chromosomal material can occur to a 
much larger extent 

Microsatellfte-based Detection of Minor Areas of Loss- 
es—In TCC 733, several chromosomal areas exhibiting 0NA 
amplification were preceded or followed by areas with a nor- 
mal CGH but reduced mRNA expression (see Fig. 1 1 TCC 733 
chromosome 1q32, 2p21, and 7q21 and q32, 9q34, and 
10q22). To determine whether these results were because of 
undetected loss of chromosomal material in these regions or 



because of other non-structural mechanisms regulating tran- 
scription, we examined two mlcrosatellites positioned at chro- 
mosome 1q25~32 and two at chromosome 2p22. Loss of 
heterozygosity (LOH) was found at both 1q25 and at 2p22 
indicating that minor deleted areas were not detected with the 
resolution of CGH (Rg. 3). Additionally, chromosome 2p in 
TCC 733 showed a CGH pattern of gairVno change/gain of 
DNA that correlated with transcript increase/decrease/in- 
crease. Thus, for the areas showing increased expression 
there was a correlation with the DNA copy number alterations 
(Fig. *\A). As indicated above, the mRNA decrease observed in 
the middle of the chromosomal gain was because of LOH, 
implying that one of the mechanisms for mRNA down-regu- 
lation may be regions that have undergone smaller losses of 
chromosomal material. However, this cannot be detected with 
the resolution of the CGH method. 

In both TCC 733 and TCC 827, the tetomeric end of chro- 
mosome 11p showed a normal ratio in the CGH analysis; 
however, clusters of five and three genes, respectively, lost 
their expression. Two mlcrosatellites (D11S1760, D11S922) 
positioned close to MUC2, IGF2, and cathepsln D Indicated 
LOH as the most likely mechanism behind the loss.of expres- 
sion (data not shown). 

A reduced expression of mRNA observed in TCC 733 at 
chromosomes 3q24, 11p11, 12p12.2, 12q21.1, and 16q24 
and in TCC 827 at chromosome 11p15.5, 12p11, 15q11.2, 
and 18q12 was also examined for chromosomal losses using 
microsateilites positioned as close as possible to the gene loci 
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Fig. 3. Microsatellrte analysis of loss of heterozygosity. Tumor 
733 showing loss of heterozygosity at chromosome 1q25, detected 
(a) by D1S215 close to Hu class I histocompatibility antigen (gene 
number 38 in Fig. 1), (b) by D1S2735 close to cathepsln E (gene 
number 41 in Fig. 1), and (c) at chromosome 2p23 by D2S2251 close 
to general 0-spectrin (gene number 1 1 on Fig. 1) and of (d) tumor 827 
showing loss of heterozygosity at chromosome 18q12 by S18S1118 
close to mitochondrial 3-oxoacyl-coenzyme A thiol ase (gene number 
12 in Fig. 1). The upper curves show the electropherogram obtained 
from normal DNA from leukocytes (A/), and the tower curves show the 
electropherogram from tumor DNA (7). In ail cases one allele is 
partially lost in the tumor ampltcon. 

showing reduced mRNA transcripts. Only the microsateliite 
positioned at 18q12 showed LOH (Fig. 3), suggesting that 
transcriptional down-regulation of genes in the other regions 
may be controlled by other mechanisms. 

Rotation between Changes in mRNA and Protein Levels- 
2D-PAGE analysis, in combination with Coomassle Brilliant 
Blue and/or silver staining, was carried out on ail four tumors 
using fresh biopsy material. 40 well resolved abundant known 
proteins migrating In areas away from the edges of the pH 
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Fig. 4. Correlation between protein levels as judged by 20- 
PAGE and transcript ratio. For comparison proteins were divided in 
three groups, unaltered in level or up- or down-regulated horizontal 
axis). The mRNA ratio as determined by oligonucleotide arrays was 
plotted for each gene {vertical axis). A, mRNAs that were scored as 
present in both tumors used for the ratio calculation; A, mRNAs that 
were scored as absent in the invasive tumors (along horizontal axis) or 
as absent in non-invasive reference (top of figure). Two different 
scalings were used to exclude scaling as a con founder, TCCs 827 
and 532 (JlA) were scaled with background suppression, and TCCs 
733 and 335 (#0) were scaled without suppression. Both compari- 
sons showed highly significant (p < 0.005) differences in mRNA ratios 
between the groups. Proteins shown were as follows: Group A (from 
Jeff), phosphoglucomutase 1 , glutathione transferase class /i number 
4, fatty acid-binding protein homologue, cytokeratin 15, and cyto- 
keratin 13; S (from te/f), fatty acid-binding protein homologue, 28-kDa 
heat shock protein, cytokeratin 13, and caleyclin; C{fromteft), «-eno- 
lase, hnRNP Bl, 28-kDa heat shock protein, 14-3-3-<; and 
pre-mRNA splicing factor; D t mesothelial keratin K7 (type II); E (from 
fop), glutathione S-transferase-Tr and mesotheflai keratin K7 (type II); 
F(from top and feft), aclenyiyi cyclase-associated protein. E-cadherin, 
keratin 19, calgizzarin. phosphoglycerate mutase, annexin IV, cy- 
toskeletal y-actln, hnRNP A1, Integral membrane protein calnexin 
(IP90). hnRNP H, brain-type dathrin light chain-a, hnRNP F, 70-kDa 
heat shock protein, heterogeneous nuclear ribonudeoprotein A/B. 
translationaliy controlled tumor protein, liver glyceraldehyde-3-phos- 
phate dehydrogenase, keratin 8, aldehyde reductase, and Na*K- 
ATPase 0-1 subunit; G, (from top and left), TCP20, calgizzarin, 70- 
kDa heat shock protein, calnexin, hnRNP H, cytokeratin 15, ATP 
synthase, keratin 19, triosephosphate teomerase, hnRNP F, liver glyc- 
eraldehyde-3-phosphatase dehydrogenase, glutathione S-transfer- 
ase-7r, and keratin 8; H (from left), plasma gelsdlin, autoantigen cal- 
retteuBn, thioredoxin. and NAD+ -dependent 15 hydroxyprostaglandin 
dehydrogenase; / (from top), prolyl 4-hydroxytase j^subunit, cyto- 
keratin 20, cytokeratin 17, prohibition, and fructose 1 ^phos- 
phatase; J annexin II; K, annexin IV; L (from top and /eft), 90-kDa heat 
shock protein, prolyl 4-hydroxyiase 0-subunit, a-enotase, GRP 78, 
cyclophilin, and cofilin. 

gradient, and having a known chromosomal location, were 
selected for analysis in the TCC pair 827/532. Proteins were 
identified by a combination of methods (see "Experimental 
Procedures"'). In general there was a highly significant corre- 
lation (p < 0.005) between mRNA and protein alterations (Fig. 
4). Only one gene showed disagreement between transcript 
alteration and protein alteration. Except for a group of cyto- 
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FfG. 5. Comparison of protein and transcript levels In invasive 
and non-invasive TCCs. The upper part of the figure shows a 2D gel 
(left) and the oligonucleotide array (rfcftf) of TCC 632. The red rectan- 
gles on the upper gel highlight the areas that are compared below. 
Identical areas of 2D gels of TCCs 532 and 827 are shown below. 
Clearly, cytokeratins 13 and 15 are strongly down-regulated In TCC 
827 (red annotation). The tile on the array containing probes for 
cytokeratin 15 is enlarged below the array (red arrow) from TCC 532 
and Is compared with TCC 827. The upper row of squares in each tife 
corresponds to perfect match probes; the lower row corresponds to 
mismatch probes containing a mutation (used for correction for un- 
specific binding). Absence of signal Is depicted as black, and the 
higher the signal the lighter the color. A high transcript level was 
detected In TCC 532 (6151 units} whereas a much lower level was 
detected in TCC 827 (absence of signals). For cytokeratin 13. a high 
transcript level was also present In TCC 532 (15659 units), and a 
much lower level was present in TCC 827 (623 units). The 2D gels at 
the bottom of the figure Qeft) show levels of PA-FABP and adipocyte- 
FABP in TCCs 335 and 733 (invasive), respectively- Both proteins are 
down-regulated in the invasive tumor. To the right we show the array 
tiles for the PA-FABP transcript A medium transcript level was de- 
tected In the case of TCC 335 (1277 units) whereas very low levels 
were detected In TCC 733 (166 units). /£F, Isoelectric focusing. 



keratins encoded by genes on chromosome 17 (Fig. 6) the 
analyzed proteins did not belong to a particular family. 26 well 
focused proteins whose genes had a know chromosomal 
location were detected in TCCs 733 and 335, and of these 19 
correlated (p < 0.005) with the mRNA changes detected using 
the arrays (Fig. 4). For example, PA-FABP was highly ex- 
pressed in the non-invasive TCC 335 but lost In the Invasive 
counterpart (TCC 733; see Fig. 5). The smaller number of 
proteins detected in both 733 and 335 was because of the 
smaller size of the biopsies that were available. 

11 chromosomal regions where CGH showed aberrations 
that corresponded to the changes in transcript levels also 
showed corresponding changes in the protein level (Table H). 
These regions included genes that encode proteins that are 
found to be frequently altered in bladder cancer, namely 
cytokeratins 17 and 20, annexins II and IV, and the fatty 
acid-binding proteins PA-FABP and FBP1. Four of these pro- 
teins were encoded by genes in chromosome 17q, a fre- 
quently amplified chromosomal area in invasive bladder 
cancers. 

DISCUSSION 

Most human cancers have abnormal DNA content, having 
lost some chromosomal parts and gained others. The present 
study provides some evidence as to the effect of these gains 
and losses on gene expression in two pairs of noo-lnvasive 
and invasive TCCs using high throughput expression arrays 
and proteomics, in combination with CGH. In general, the 
results showed that there is a clear individual regulation of the 
mRNA expression of single genes, which in some cases was 
superimposed by a DNA copy number effect In most cases, 
genes located in chromosomal areas with gains often exhib- 
ited increased mRNA expression, whereas areas showing 
losses showed either no change or a reduced mRNA expres- 
sion. The latter might be because of the fact that losses most 
often are restricted to loss of one allele, and the cut-off point 
for detection of expression alterations was a 2-fold change, 
thus being at the border of detection. In several cases, how- 
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* In cases where the corresponding alterations were found in both TCCs 827 and 733 these are shown as 827/733. 



Molecular & Cellular Proteomics 1. 1 43 



Gene Copy Numbers, Transcripts, and Protein Levels 



ever, an increase or decrease in DNA copy number was 
associated with de novo occurrence or complete loss of tran- 
script, respectively. Some of these transcripts could not be 
. detected in the non-invasive tumor but were present at rela- 
tively high levels in areas with DNA amplifications in the inva- 
sive tumors (e.g. in TCC 733 transcript from cellular iigand of 
annexin II gene (chromosome 1q21) from absent to 2670 
arbitrary ttnits; in TCC 827 transcript from small proline-rich 
protein 1 gene (chromosome 1q12-q21.1) from absent to 
1326 arbitrary units). It may be anticipated from these data 
that significant clustering of genes with an increased expres- 
sion to a certain chromosomal area indicates an increased 
likelihood of gain of chromosomal material in this area 

Considering the many possible regulatory mechanisms act- 
ing at the level of transcription, it seems striking that the gene 
dose effects were so clearly detectable in gained areas. One 
hypothetical explanation may tie In the loss of controlled 
methylation In tumor cells (17-19). Thus, it may be possible 
that in chromosomes with increased DNA copy numbers two 
or more alleles coutd be demethylated simultaneously leading 
to a higher transcription level, whereas in chromosomes with 
losses the remaining allele could be partly methylated, turning 
off the process (20, 21). A recent report has documented a 
pioidy regulation of gene expression in yeast, but in this case all 
the genes were present in the same ratio (22), a situation that is 
not analogous to that of cancer ceils, which show marked 
chromosomal aberrations, as well as gene dosage effects. 

Several CGH studies of bladder cancer have shown that 
some chromosomal aberrations are common at certain 
stages of disease progression, often occurring in more than 1 
of 3 tumors. In pTa tumors, these include 9p- f 9q 1 q + , Y- 
(2, 6), and UvpTI tumors, 2q-,11p-. 11q-, 1q+, 5p+, 8q+, 
17q+, and 20q+ (2-4, 6, 7). The pTa tumors studied here 
showed similar aberrations such as 9p- and 9q22-q33- and 
9q- and Y-, respectively. Likewise, the two minimal invasive 
pT1 tumors showed aberrations that are commonly seen at 
that stage, and TCC 627 had a remarkable resemblance to the 
commonly seen pattern of losses and gains, such as 1 q22-24 
amplification (seen in both tumors), 1 1q14-q22 loss, the latter 
often linked to 17 q+ (both tumors), and 1q+ and 9p~, often 
linked to 20q+ and 11 q13+ (both tumors) (7-9). These ob- 
servations Indicate that the pairs of tumors used In this study 
exhibit chromosomal changes observed in many tumors, and 
therefore the findings could be of general. Importance for 
bladder cancer 

Considering that the mapping resolution of CGH is of about 
20 megabases it is only possible to get a crude picture of 
chromosomal Instability using this technique* Occasionally, 
we observed reduced transcript levels close to or inside re- 
gions with increased copy numbers. Analysis of these regions 
by positioning heterozygous microsateliites as close as pos- 
sible to the locus showing reduced gene expression revealed 
loss of heterozygosity in several cases. It seems likely that 
multiple and different events occur along each chromosomal 



arm and that the use of cDNA microarrays for analysis of DNA 
copy number changes will reach a resolution that can resolve 
these changes, as has recently been proposed (2). The outlier 
data were not more frequent at the boundaries of the CGH 
aberrations. At present we do not know the mechanism be- 
hind chromosomal aneuploldy and cannot predict whether 
chromosomal gains will be transcribed to a larger extent than 
the two native alleles. A mechanism as genetic imprinting has 
an impact on the expression level In normal cells and is often 
reduced in tumors. However, the relation between imprinting 
and gain of chromosomal material is not known. 

We regard it as a strength of this investigation that we were 
able to compare invasive tumors to benign tumors rather than 
to normal urothelium, as the tumors studied were biologically 
very close, and probably may represent successive steps in 
the progression of bladder cancer. Despite the limited amount 
of fresh tissue available it was possible to apply three different 
state of the art methods. The observed correlation between 
DNA copy number and mRNA expression is remarkable when 
one considers that different pieces of the tumor biopsies were 
used for the different sets of experiments. This indicate that 
bladder tumors are relatively homogenous, a notion recently 
supported by CGH and LOH data that showed a remarkable 
similarity even between tumors and distant metastasis (10, 23). 

In the few cases analyzed, mRNA. and protein levels 
showed a striking correspondence although in some cases 
we found discrepancies that may be attributed to translationai 
regulation, post-translatlonal processing, protein degrada- 
tion, or a combination of these. Some transcripts belong to 
under-translated mRNA pools, which are associated with few 
transiationaily inactive ribosomes; these pools, however, 
seem to be rare (24). Protein degradation, for example, may 
be very Important in the case of polypeptides with a short 
half-life (e.g. signaling proteins). A poor correlation between 
mRNA and protein levels was found in liver cells as deter- 
mined by arrays and 2D-PAGE (25), and a moderate correla- 
tion was recently reported by Ideker et al. (26) in yeast 
(Interestingly, our study revealed a much better correlation 
between gained chromosomal areas and Increased mRNA 
levels than between toss of chromosomal areas and reduced 
mRNA levels. In general, the level of CGH change determined 
the ability to detect a change in transcript?) One possible 
explanation could be that by losing one allele the change in 
mRNA level is not so dramatic as compared with gain of 
material, which can be rather unlimited and may lead to a 
severatfold Increase in gene copy number resulting In a much 
higher Impact on transcript level, the latter would be much 
easier to detect on the expression arrays as the cut-off point 
was placed at a 2 -fold level so as not to be biased by noise on 
the array. Construction of arrays with a better signal to noise 
ratio may In the future allow detection of lesser than 2-fold 
alterations in transcript levels, a feature that may facilitate the 
analysis of the effect of loss of chromosomal areas on tran- 
script levels. 
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In eleven cases we found a significant correlation between 
DNA copy number, mRNA expression, and protein level. Four 
of these proteins were encoded by genes located at a fre- 
quently amplified area in chromosome 17q. Whether DNA 
copy number Is one of the mechanisms behind alteration of 
these eleven proteins is at present unknown and will have to 
be proved by other methods using a larger number of sam- 
ples. One factor making such studies complicated is the targe 
extent of protein modification that occurs after translation, 
requiring immunoidentiflcation and/or mass spectrometry to 
correctly Identify the proteins in the gels. 

In conclusion, the results presented in this study exemplify 
the large body of knowledge that may be possible to gather in 
the future by combining state of the art techniques that follow 
the pathway from DNA to protein (26). Here, we used a tradi- 
tional chromosomal CGH method, but in the future high reso- 
lution CGH based on microanrays with many thousand radiation 
hybrid-mapped genes will increase the resolution and informa- 
tion derived from these types of experiments (2). Combined with 
expression arrays analyzing transcripts derived from genes with 
known locations, and 2D gel analysis to obtain information at 
the post-translational level, a clearer and more developed un- 
derstanding of the tumor genome will be forthcoming. 
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ABSTRACT 

Genetic changes underlie -tumor progression and may lead to cancer- 
specific expression of critical genes. Over 110Q publications have de- 
; scribed the use of comparative genomic hybridization (CGH) to analyze 
the pattern of copy number alterations in cancer, but very few of the genes 
affected are known. Here, we performed high-resolution CGH analysis on 
cDNA micro arrays In breast cancer and directly compared copy number 
and mRNA expression levels of 13,824 genes to quantifete the impact of k 
genomic changes on gene expression. We identified and mapped the 
boundaries of 24 independent amplicons, ranging In size from 0 J to 12 
Mb, Throughout the genome, both high- and low-level copy number 
changes had a substantial impact on gene expression* with 44% of the 
highly amplified genes showing overexpression and 10.5% of me highly 
overexpressed genes being amplified. Statistical analysis with random 
permutation tests identified 270 genes -whose expression levels across' 14 
samples were systematically attributable to gene, amplification. TJiese 
included most previously described amplified genes in breast cancer and 
many novel targets for genomic alterations, including the BOXB7 gene, 
the presence of which in a novel amplicon at 17q21J was validated in 
10.2% of primary breast cancers and associated with poor patient prog- 
nosis. In conclusion,' CGH on cDNA microarrays revealed hundreds of • 
novel genes whose overexpression Is attributable to gene amplification. 
These genes may provide insights to the clonal evolution and progression 
of breast cancer and highlight promising therapeutic targets. 

INTRODUCTION 

Gene expression patterns revealed by cDNA microarrays have 
facilitated classification of cancers into biologically distinct catego- 
ries, some of which may explain the clinical behavior of the tumors 
(1-6). Despite this progress in diagnostic classification, the molecular 
mechanisms underlying gene expression patterns in cancer have re- 
rnainedi! elusive, and the utility of gene expression profiling in the 
identification of specific tnerapeufic targets remains ImiSSdr^ 

Accumulation of genetic defects is thought to underlie the clonal 
evolution of cancer. Identification of the genes that mediate the effects 
of genetic changes may be important by highlighting transcripts that 
are actively involved in tumor progression. Such transcripts and their 
encoded proteins would be ideal targets , for anticancer therapies, as 
demonstrated by the clinical success of new therapies against ampli- 
fied oncogenes, such as ERBB2 md EGFR (J 9 8), in breast cancer and 
other solid tumors. Besides amplifications of known oncogenes, over 
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Fig. 1. Impact of gene copy number on global gene expression levels. A, percentage of 
over- and uno^ata^pressedqgenes (Y axis) according to copy .number^atios (<X axis). 
Threshold values used for over- and underexpressfon were >2J84 (global upper 7% of 
the cDNA ratios) and <0.4826 (global lower 7% of the expression ratios). B, percentage 
of amplified and deleted genes according to expression ratios. Threshold values for 
amplification and deletion were > 1 .5 and <0.7. *' 



20. recurrent regions of DNA amplification have been mapped in 
breast cancer by CGH 3 (9, 10)., However, these amplicons are often 
large and poorly defined, and their impact on gene expression remains 
unknown 

We hypothesized that genome-wide identification of those gene 
expression changes that are attributable to underlying gene copy 
number alterations would highlight transcripts that are. actively in- 
volved in the causation or maintenance of the malignant phenotype. 
To identify such transcripts, we applied a combination of cDNA and 
CGH microarrays to: (a) determine the global impact that gene copy 
number variation plays in breast cancer developmentand progression; 
and (b) identify and characterize those genes whose mRNA expres- 



3 The abbreviations used are: CGH, comparative genomic hybridization; FISH, fluo- 
rescence in situ hybridization; RT-PCR, reverse transcription-PCR. 
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sion is most significantly associated with amplification of the corre- 
sponding genomic template. 

MATERIALS AND METHODS 

Breast Cancer Oifl Lines* Fourteen breast cancer cell lines' (BT-20, BT- 
474, HOC1428, Hs578t, MCF7, KfDA-361, MDAt43.6, MDA-453 t MDA-468, 
SKBFt-3, T-47D, UACC812, ZR-75-1, and ZR-75-30) were obtained from the 
American Type Culture Collection (Manassas, VA). Cells were grown under 
recommended culture conditions. Genomic DNA and raRNA were isolated 
using standard protocols. 

Copy Number and Expression Analyses by cDNA Microarrays. The 
preparation and printing of the 13,824 cDNA clones on glass slides were 
performed as described (1 1-13). Of these clones, 244 represented uncharac- 
terized expressed sequence tags, and the remainder corresponded to known * 
genes. CGH experiments on cDNA rru'eroarrays were done as described (14, 
15), Briefly, 20 jtg of genomic DNA from breast cancer cell lines and normal 
human WBCs were digested for 14.-48 h with Ahjl-wd Asal (Life Technol- 
ogies, Inc., Rockville, MD) and purified by phenol/chloroform extraction. Six 
Mg of digested cell line DNAs were labeled with Cy3-dUTP (Amersham 
Pharmacia) and normal DNA with Cy5*dUTP (Amersham Pharmacia) using 
the Bioprime Labeling kit (Life Technologies, Inc.). Hybridization (14, 15) and. 
posthybridization washes (13) were done as. described For the expression 
analyses, a standard reference (Universal Human Reference RNA; Stratagene, 
La Jolla, CA) was used in all experiments. Forty jtg of reference RNA were 
labeled with Cy3-dUTP and 3.5 jxg of test. mRNA with Cy5.-d!JTP, and the 
labeled cDNAs were hybridized on microarrays as described (13, 15). For both 
microarray analyses, a laser confocal scanner (Agilent .Technologies, Palo 
Alto, CA) was used to measure the fluorescence intensities at the target 
locations using the 0EARRAY software (16). After background subtraction, 
average intensities at each clone in the test hybridization were divided by the 
average intensity of me corresponding clone in the control hybridization. For 
the copy number analysis, the ratios were normalized on the basis of the 
distribution of ratios of all targets on the array and for the expression analysis 
on the basis of 88 housekeeping genes, which were spotted four times onto the 
array. Low quality measurements (/.*, copy number data with mean reference 
intensity <100 fluorescent units, and expression data with both test , and 
reference intensity <100 fluorescent units and/or with spot size <50 units) 



were excluded from the analysis and were treated as missing values. The 
distributions of fluorescence ratios were used to define cutpoints for increased/ 
decreased copy number. Genes with CGH ratio >M3 (representing the upper 
5% of the CGH ratios across all experiments) were considered to be amplified, 
and genes with ratio <0:73 (representing the lower 5%) were considered to be 
deleted. 

Statistical Analysis of CGH and cDNA Microarray Data. To evaluate 
the influence of copy number alterations on gene expression, we. applied the 
following statistical approach. CGH and cDNA calibrated intensity ratios were 
log-transformed and normalized using median centering of the values in each 
cell line. Furthermore, cDNA ratios for each gene across all 14 cell lines were 
median centered. For each gene, the CGH data were represented by a vector 
that was labeled 1 for amplification (ratio, >1.43) and 0 for no amplification. 
Amplification was correlated with gene expression using the signal-to-noise 
^statistics (1). We calculated a weight, w r for each gene as^llows: 

* 8l +.cr f0 

where m gU <r gl and denote the means and SDs for the expression 

levels for amplified and nonamplified cell lines, respectiveiy. To assess the 
statistical significance of each weight, we performed 10,000 random permu- 
tations of the label vector. The probability mat a gene had a larger or equal 
weight by random permutation than the original weight was denoted by a. A 
low a (<0,05) indicates a strong association between gene expression and 
amplification. 

Genomic Localization of cDNA Clones and Amplicon Mapping, Each 
cDNA clone on the microarray was assigned to a Unigene cluster using the 
Unigerie Build 14 J * A database of genomic sequence alignment information 
for mRNA sequences was created from die August 2001 freeze of the Uni- 
versity of California Santa Cruz's GoldenPath database. 7 The chromosome and 
bp positions for each cDNA clone were men retrieved by relating these data 
sets, Amplicons were defined as a CGH copy number ratio >2.0 in. at least two 
adjacent clones in two or more cell lines or a CGH ratio >2.0 in at least three 
adjacent clones in a single cell line. The amplicon start and end positions were 



* Internet address: http://research.iihgri.nik^ cdna.htmL 
7 Internet address: www.genomc.ucsc.edu. ~ 
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Table 1 Summary of independent amplicons in J4 breast cancer cell lines by 

CGH microarray 

Location 

lpI3 
lq21 
lq22 
3pI4 

7pl2.1-7pll,2 
7q31 
7q32 

8q21.11-8q2U3 
8q213 

8q233-*q24.14 
8q24.22 
9pl3 

13q22-q31 . 
I6a22 
17qll 

17qlZ-q2l.2 
17o21.32-q2l.33 
17q22-q23.3 
17o233-<£4.3 
19ql3 
20qll.22 
20ql3.12 
20ql3.12-<il3,13 
20ql3^-ql332 



Start (Mb) 


End (Mb) 


Size (Mb) 


|M 70 


132.94 


.0.2 


ITt <» 

1 ilCfL 


If*! ft 

177.25 


33 


170 70 


179.57 


03 


'71 OA 


74.66 . 


2.7 




60.95 


53 


1 23.73 


130*96 


5.2 


1W.U1 


140.68 


0.7 


86.45 ' 


92.46 


6.0 


98.45 


103.05 


4.6 


129.88 


142.15 


123 


.15121 


152.16 


1.0 


38.65 


.39.25. 


0.6 


77.15 


81.38 


4.2 


86.70 


87.62 


0.9 


29.30 


30L85 


1.6 


39.79 


42.80 


3.0 


52.47 


55.80 


33 


63.81 


69.70 


5.9 


69.93 


74.99 


5.1 


40.63 


41.40 


0.8 


.34.59 


35.85 


.13 


44.00 


45.62 


1.6 


46.45 


49.43 


3.0 


502 


59.12 


. . 7.8 



CGH were validated, with lq21, 17ql2-q21.2* 17q22^q23, 20ql3.l 
and 20ql3.2 regions being most commonly amplified. Furthermore' 
the boundaries of these amplicons were precisely delineated. In ad - 
dition, novel amplicons were identified at 9pl3 (38.65-39 25 MbV 
and 17q2L3 (52.47-55,80 Mb). " • 

Direct Identification of Putative Amplification Target Genes. 
The cDNA/CGH microarray technique enables the direct correla- 
tion of copy number and expression data on a gene-by-gene basis 
throughout the genome. We directly annotated hi^resolution 
CGH plots with gene expression data using color coding. Fig. 2C 
shows that most of the amplified genes hi the MCF-7 breast cancer 
cell line at lpl3, 17q22-q23, and 20ql3 were highly overex- 
pressed. A view of chromosome 7 in the- MDA-468 cell line 
implicates EGFR as the most highly overexpressed and amplified 
gene at 7pl l-p!2 (Fig. 3,4). In BT-474, the two known amplicons 
at 17ql2 and 17q22-q23 contained numerous highly overex- 
pressed genes (Fig. 3B). In addition, several genes, including the 
homeobox genes HOXB2 and HOXB7, were highly amplified in a 
. previously undescribed independent amplicon. at. 17q2 1.3. HOXB7 
was systematically amplified (as validated by FISH, Fig. %B, insef) 
as well as overexpressed (as verified:by RT-PCR, data not shown) 
in BT-474, UACC812, and ZR-75-30 cells. Furthermore, this novel 



extended to include neighboring nonamplified clones (ratio, <15). The am- 
plicon size determination was partially dependent on local clone density. 

FISfiL Dual-color interphase FISH to breast cancer cell lines Was done as 
described (17). Bacterial artificial chromosome clone RP11-361K8 was la- 
beled with Sr>ectnimOrange (Vysis, Downers Grove, IL), and Spectrum- 
Orange-labeled probe for EGFR was obtained from Vysis. SpectrumGreen- 
labeled chromosome 7 and 17 centromere probes (Vysis) were used as a 
reference. A tissue microarray containing 612 fonnalin-fixed, paraffin-embed- 
ded primary breast cancers (17) was applied in FISH analyses as described 
, (18). The use of these specimens was approved by the Ethics Committee of the 
University of Basel and by the NIH. Specimens containing a 2-fold or higher 
increase in the number oftest probe signals, as compared with cwresponding 
centromere, signals, in at least 10% of the tumor cells were considered to be 
amplified. Survival analysis was performed using the Kaplan-Meier method 
and the log-rank test 

RT-PCR. The HOXB7 expression level was dWmined relative to 
GAPDH. Reverse transcription and PGR amplification were performed using 
Access RT-PCR System (Promega Corp., Madison, WI) with 10 ng ofmRNA 
as a template. HOXB7 primers were 5 '^AGCAGAGGGACTCGGACTT-3 ' 
and 5'-GCGTCAGOTAGCG ATTOTAG-3 ' . 

RESJJLTS 

Global Effect of Copy Number on Gene Expression. 13,824 
arrayed cDNA clones were applied for analysis of gene expression 
and gene copy number (CGH micfoarrays) in 14 breast cancer cell 
lines. The results illustrate a considerable influence of copy number 
on gene expression patterns. Up to 44% of the highly amplified 
transcripts (CGH ratio, >2.5> were overexpressed (Le., belonged to 
the global upper 7% of expression ratios), compared with only 6% for 
genes with normal copy number levels (Fig. 1 A). Conversely, 10.5% 
of the transcripts, with high-level expression (cDNA ratio, >10) 
showed increased copy number (Fig. IB). Low-level copy number 
increases and decreases were also associated with similar, although 
less dramatic, outcomes on gene expression (Fig. 1). 

Identification of Distinct Breast Cancer Amplicons* Base-pair 
locations obtained for 1 1,994 cDNAs (86.8%) were used to plot copy 
number changes as a function of genomic position (Fig. 2, Supple- 
ment Fig. A). The average spacing of clones throughout the genome 
was 267 kb; This high-resolution mapping identified 24 independent 
breast cancer amplicons, spanning from 02 to-12 Mb of DNA (Table 
1). Several amplification sites detected previously by chromosomal 





Fig. 3. Annotation of gene expression data on CGH microarray profiles. A, genes in the 
Tpll -pl2 amplicon in the MDA-468 cell line are highly expressed {red dots) and include 
the EGFR oncogene, B, several genes in Ac I7qi2, . 17q2U, and 17q23 amplicons in the 
BT-474 breast cancer cell line are highly overexpressed (red) and include the HOXB7 
gene. The data labels and color coding are as indicated for Fig. 2C Insets show 
chromosomal CGH profiles for the corresponding chromosomes and validation of the 
increased copy number by interphase FISH using EGFR (red) and chromosome 7 
centromere probe (green) to MDA-468. (4) and HOX87*pcc\fic probe (red) and chro- 
mosome 17 centromere (green) to BT-474 ceils (B^ 
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level copy number increase. Low-level copy number gains and losses 
also had a significant influence on expression levels of genes in the 
regions affected, but these effects were more subtle on a gene-by-gene 
basis than those of high-level amplifications. However, the impact of 
low-level gains on the dysregulation of gene expression patterns in 
cancer may be equally important if not more important than that of 
high-level amplifications. Aneuploidy and low-level gains and losses 
of chromosomal arms represent the most common types of genetic 
alterations in breast and other cancers and, therefore, have an influ- 
ence on many genes. Our results in breast cancer extend the recent 
studies on the impact of aneuploidy on global gene expression pat-" 

.. terns in yeast cells, acute myeloid leukemia, and a prostate cancer 
model system (22-24). 

The CGH microarray analysis identified 24 independent breast 
cancer amplipons. We defined the precise boundaries for, many am- 

. plfcons detected previously by chromosomal CGH (?, 10, 25, 26) and 

, also discovered novel amplicons that had not been detected previ- 
ously, presumably because of their small size (only 1-2 Mb) or close 
proximity to other larger amplicons. One of these novel amplicons 
involved the homeobox gene region at 17q21 3 and led to the over- 
expression of the HOXB7 and HOXB2 genes. Hie homeodomain 
transcription factors are"' known to be key regulators of embryonic 
development and have been occasionally reported to undergo aberrant 
expression in cancer (27, 28). HOXB7 transaction induced cell pro- 
liferation in melanoma, breast, and ovarian cancer cells and increased 
tumorigenicity and angiogenesis in breast cancer (29-32), the pres- 
ent results imply that gene amplification may be a prominent mech- 
anism for overexpressing HOXB7 in breast cancer and suggest that 
HOXB7: contributes to tumor progression and confers an aggressive 
disease phenotype in breast cancer. This view is supported by our 
finding of amplification oi HOXB7 in 10% of 363 primary breast 
cancers, as well as an association of amplification with poor prognosis 

: of the patients. 

We carried out a systematic search to identity genes whose 
expression levels across all 14 cell lines were attributable to 
amplification status.. Statistical analysis revealed 270 such genes 
(representing -2% of all genes pn the axray), including not only 
previously . described amplified genes, such as HER-2, MYC t 
EGFR, ribosomal protein s6 kinase, and AIB3, but also numerous 
novel genes such as NRAS*reIated gene (lpl3), syndecan-2 (8q22), 
and bone morphogentc protein (20ql3.1), whose activation by 
amplitlcatiqp may similarly promote breast cancer progression. 
Most of the 270. genet have not been implicated previously in 
breast cancer development and suggest novel pathogenetic mech- 
anisms. Although we would not expect all of them to be causally 
involved, it is intriguing that 84% of the genes with associated 
functional information were implicated in apoptosis, cell prolifer- 
ation, signal transduction, transcription, or other cellular processes 
that could directly imply a possible role in cancer progression. 
Therefore, a detailed characterization of these genes may provide 
biological insights to breast cancer progression and might lead to 
the development of novel therapeutic strategies. 

In summary, we. demonstrate application of cDNA microarrays 
to the analysis of both copy number and expression levels of over 
12,000 transcripts throughout the breast cancer genome, roughly 
once every 267 kb. This analysis provided: (a) evidence of a 
,j, prominent global influence of copy number changes on gene 
expression levels; (b) a high-resolution map. of 24 independent 
amplicons in breast cancer, and (c) identification of a set of 270 
genes, the overexpression of which was statistically attributable to 
gene amplification; Characterization of a novel amplicon at 
17q21.3 implicated amplification and overexpression of the 
HOXB7 gene in breast cancer, including a clinical association 
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between HOXB7 amplification and poor patient prognosis. Overall 
our results illustrate how the identification of genes activated by 
gene amplification provides a powerful approach to highlight 
genes with an important role in cancer as well as to prioritize and 
validate putative targets for therapy development 
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Microarray analysis reveals a major direct role of 
DNA copy number alteration in the transcriptional 
program of human breast tumors 
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Genomic DNA copy number alterations are key genetic events in 
the development and progression of human cancers. Here we 
report a genome-wide microarray comparative genomic hybrid- 
feation (array CGH) analysis of DNA copy number variation in 
a series of primary human breast tumors. We have profiled DNA 
copy number alteration across 6,691 mapped human genes, in 44 
predominantly advanced, primary breast tumors and 10 breast 
cancer ceti lines. White the overall patterns of DNA amplification 
and deletion corroborate previous cytogenetic studies, the high- 
resolution (gene-by-gene) mapping of amplicon boundaries and 
the quantitative analysis of amplicon shape provide significant 
improvement In the localization of candidate oncogenes. Parallel 
microarray measurements of mRNA levels reveal the remarkable 
degree to which variation in gene copy number contributes to 
variation in gene expression in tumor cells. Specifically, we find 
that 62% of highly amplified genes show moderately or highly 
elevated expression, that DNA copy number influences gene ex- 
pression across a wide range of DNA copy number alterations 
(deletion, low-, mid- and high-level amplification), that on average, 
a 2-fold change in DNA copy number Is associated with a corre- 
sponding 1.5-fold change in mRNA levels, and that overall, at least 
12% of all the variation in gene expression among the breast 
tumors Is directly attributable to underlying variation in gene copy 
number. These findings provide evidence that widespread DNA 
copy number alteration can lead directly to global deregulation of 
gene expression, which may contribute to the development or 
progression of cancer, 

Conventional cytogenetic techniques, including comparative 
genomic hybridization (CGH) (1), have led to the identifi- 
cation of a number of recurrent regions of DNA copy number 
alteration in breast cancer cell lines and tumors (2-4). While 
some of these regions contain known or candidate oncogenes 
[e.g., FOFR1 (8pll), MYC (8q24), CCND1 (Uql3), ERBB2 
(r7ql2), and ZNF217 (20ql3)] and tumor suppressor genes 
(RBI (13ql4) and TP53 (Hptf)], the relevant gene(s) within 
other regions (e.g, ( gain of lq, 8q22, and 17q22-24, and loss of 
8p) remain to be identified. A high-resolution genome-wide 
map, delineating the boundaries of DNA copy number alter- 
ations in tumors, should facilitate the localization and identifi- 
cation of oncogenes and tumor suppressor genes in breast 
cancer. In this study, we have created such a map, using 
array-based CGH (5-7) to profile DNA copy number alteration 
in a series of breast cancer cell lines and primary tumors. 

An unresolved question is the extent to which the widespread 
DNA copy number changes that we and others have identified 
in breast tumors alter expression of genes within involved 
regions. Because we had measured mRNA levels in parallel in 
the same samples (8) t using the same DNA microarrays, we had 
an opportunity to explore on a genomic scale the relationship 
between DNA copy number changes and gene expression. From 



this analysis, we have identified a significant impact of wide- 
spread DNA copy number alteration on the transcriptional 
programs of breast tumors. 

Materials and Methods 

Tumors and Cell Lines, Primary breast tumors were predominantly 
large (>3 cm), intermccbate-grade, infiltrating ductal carcino- 
mas, with more than 50% being lymph node positive. The 
fraction of tumor cells within specimens averaged at least 50%. 
Details of individual tumors have been published (8, 9), and 
are summarized in Table 1, which is published as supporting 
information on the PNAS web site, www.pnas.org. Breast cancer 
cell lines were obtained from the American Type Culture 
Collection. Genomic DNA was isolated either using Qiagen 
genomic DNA columns, or by phenol/chloroform extraction 
followed by ethanol precipitation. 

DNA Labeling and Microarray Hybridizations. Genomic DNA label- 
ing and hybridizations were performed essentially as described 
in Pollack et al (7), with slight modifications. Two micrograms 
of DNA was labeled in a total volume of 50 microliters and the 
volumes of all reagents were adjusted accordingry, "Test" DNA 
(from tumors and cell lines) was f luorescentiy labeled (Cy5) and 
hybridized to a human cDNA microarray containing 6,691 
different mapped human genes (i.e., UniGene clusters). The 
reference" (labeled with Cy3) for each hybridization was nor- 
mal female leukocyte DNA from a single donor. The fabrication 
of cDNA microarrays and the labeling and hybridization of 
mRNA samples have been described (8). 

Data Analysis and Map Positions. Hybridized arrays were scanned 
on a GenePix scanner (Axon Instruments, Foster City, CA), and 
fluorescence ratios (test/reference) calculated using scan alyze 
software (available at http://rana.lbl.gov). Fluorescence ratios 
were normalized for each array by setting the average log 
fluorescence ratio for all array elements equal to 0. Measure- 
ments with fluorescence intensities more than 20% above back- 
ground were considered reliable. DNA copy number profiles 
that deviated significantly from background ratios measured in 
normal genomic DNA control hybridizations were interpreted as 
evidence of real DNA copy number alteration (see Estimating 
Significance of Altered Fluorescence Ratios in the supporting i 
information). When indicated, DNA copy number profiles are 
displayed as a moving average (symmetric 5-nearest neighbors), I 
Map positions, for arrayed human cDNAs were assigned by 
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identifying the starting position of the best and longest match of 
any DNA sequence represented in the corresponding UniGene 

^ C // (10) against the MGold « n Path" genome assembly 
(http://genome.ucsc.edu/; Oct 7, 2000 Freeze). For UniGene 
clusters represented by multiple arrayed elements, mean fluo- 
rescence ratios (for all elements representing the same UniGene 
cluster) are reported. For mRNA measurements, fluorescence 
ratios are "mean-centered" (Le„ reported relative to the mean 
ratio across the 44 tumor samples). The data set described here 
can be accessed in its entirety in the supporting information. 

Results 

We performed CGH on 44 predominantly locally advanced 
P^^y breast tumors and 10 breast cancer cell lines, using 
cDNA microarrays containing 6,691 different mapped human 
genes (Fig. la; also see Materials and Methods for details of 
microarray hybridizations). To take full advantage of the im- 
proved spatial resolution of array CGH, we ordered (f luores- 
cence ratios for) the 6,691 cDNAs according to the "Golden 
Path" (http://genome.ucsc.edu/) genome assembly of the draft 
human genome sequences (11). In so doing, arrayed cDNAs not 
only themselves represent genes of potential interest (e.g. 
candidate oncogenes within amphcons), but also provide precise 
genetic landmarks for chromosomal regions of amplification and 

12964 | www-pna$.org/cgi/dol/tai073/pnas. 162471999 



deletion. Parallel analysis of DNA from cell lines containing 
different numbers of X chromosomes (Fig. lb), as we did before 
(7), demonstrated the sensitivity of our method to detect single* 
copy loss (45, XO), and IS- (47 f XXX), % (48JCXXX), or 
2,5-fold (49,XXXXX) gains (also see Fig. 5, which is published 
as supporting information on the PNAS web site). Fluorescence 
ratios were linearly proportional to copy number ratios, which 
were slightly underestimated, in agreement with previous ob- 
servations (7). Numerous DNA copy number alterations were 
evident in both the breast cancer cell lines and primary tumors 
(Fig. la), detected in the tumors despite the presence of euploid 
non-tumor cell types; the magnitudes of the observed changes 
were generally lower in the tumor samples. DNA copy-number 
alterations were found in every cancer cell line and tumor, and 
on every human chromosome in at least one sample. Recurrent 
regions of DNA copy number gain and loss were readily iden- 
tifiable. For example, gains within Iq, 8q, 17q, and 20q were 
observed in a high proportion of breast cancer cell lines/tumors 
(9096/69%, 100%/47%, 100%/60%, and 90%/44%, respective* 
ry), as were losses within lp, 3p, 8p, and 13q (80%/24% 
80%/22%, 80%/22%, and 70%/18%, respectively), consistent 
with published cytogenetic studies (refe. 2-4; a complete listing 
of gams/losses is provided in Tables 2 and 3, which are published 
as supporting information on the PNAS web site). The total 
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Rg.2. DNAcopynumberarteratlonaaossch^ 

of X chromosomes, for breast cancer cell lines, and for breast tumors. Breast cancer cell Bnes and tumors are separately ordered by hierarchical clustering to 
highlight recurrent copy number changes. The 241 genes present on the mlcroanray* and mapping to chromosome 8 are ordered by position along the 
chromosome. Fluorescence ratios (test/reference) are depicted by a log, pseudocolor scale (indicated). Selected genes are Indicated with color-coded text (red 
Increased; green, decreased; black, no change; gray, not well measured) to reflect correspondingly altered mRNA levels (observed In the majority of the subset 
of samples displaying the DNA copy number change). The map positions for genes of interest that are not represented on the microarray are Indicated in the 
row above those genes represented on the array. <b) Graphical display of DNA copy number profile for breast cancer cell line SKBR3. Ruorescence ratios 
(tumor/normal) are plotted on a log 2 scale for chromosome 8 gene*, ordered along the chromosome. 



number of genomic alterations (gains and losses) was found to 
be significantly higher in breast tumors that were high grade (P » 
0.008), consistent with published CGH data (3), estrogen recep- 
tor negative (P » 0,04), and harboring TP53 mutations (P - 
0.0006) (see Table 4, which is published as supporting informa- 
tion on the PNAS web site). 

The improved spatial resolution of our array CGH analysis is 
illustrated for chromosome 8, which displayed extensive DNA 
copy number alteration in our series. A detailed view of the 
variation in the copy number of 241 genes mapping to chromo- 
some 8 revealed multiple regions of recurrent amplification; 
each of these potentially harbors a different known or previously 
uncharacterized oncogene (Fig. 2a), The complexity of amplicdn 
structure is most easily appreciated in the breast cancer cell line 
SKBR3. Although a conventional CGH analysis of 8q in SKBR3 
identified only two distinct regions of amplification (12), we 
observed three distinct regions of high-level amplification (la- 
beled 1-3 in Fig. 2b). For each of these regions we can define the 



boundaries of the interval recurrently amplified in the tumors we 
examined; in each case, known or plausible candidate oncogenes 
can be identified (a description of these regions, as well as the 
recurrently amplified regions on chromosomes 17 and 20, can be 
found in Figs. 6 and 7, which are published as supporting 
information on the PNAS web site). 

For a subset of breast cancer cell lines and tumors (4 and 37* 
respectively), and a subset of arrayed genes (6,095), mRNA 
levels were quantitatively measured in parallel oy using cDNA 
microarrays (8). The parallel assessment of mRNA levels is 
useful in the interpretation of DNA copy number changes. For 
example, the highly amplified genes that are also highly ex- 
pressed are the strongest candidate oncogenes within an ampli- 
con. Perhaps more significantly, our parallel analysis of DNA 
copy number changes and mRNA levels provides us the oppor- 
tunity to assess the global impact of widespread DNA copy 
number alteration on gene expression in tumor cells. 

A strong influence of DNA copy number on gene expression 
is evident in an examination of the pseudocolor representations 
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Rg. 3. Concordance between DNA copy number and gene expression across chromosome 1 7. DNA copy number alteration (Upper) and mRNA levels (lower) 
are illustrated lor breast cancer cell lines and tumors. Breast cancer cell lines and tumors are separately ordered by hierarchical clustering (Upper), and the 
Identical sample order is maintained Clowe/). The 354 genes present on the mkroarrays and mapping to chromosome 1 7, and for which both DNA copy number 
and mRNA levels were determined, are ordered by position along the chromosome; selected genes are Indicated in color-coded text (see Rg. 2 legend). 
Fluorescence ratios (test/reference) are depicted by separate iog 2 pseudocolor scales (indicated). 



of DNA copy number and mRNA levels for genes on chromo- 
some 17 (Fig. 3), The overall patterns of gene amplification and 
elevated gene expression are quite concordant; Le., a significant 
fraction of highly amplified genes appear to be correspondingly 
highly expressed. The concordance between high-level amplifi- 
cation and increased gene expression is not restricted to chro- 
mosome 17* Genotmvwide, of 117 high-level DNA amplifica- 
tions (fluorescence ratios >4, and representing 91 different 
genes)* 62% (representing 54 different genes; see Table 5, which 
is published as supporting information on the PNAS web site) 
are found associated with at least moderately elevated mRNA 
levels (mean-centered fluorescence ratios >2), and 42% (rep- 
resenting 36 different genes) are found associated with compa- 
rably highly elevated mRNA levels (mean-centered fluorescence 
ratios >4), 

To determine the extent to which DNA deletion and lower- 
level amplification (in addition to high-level amplification) are 
also associated with corresponding alterations in mRNA levels, 
we performed three separate analyses on the complete data set 
(4 cell lines and 37 tumors, across 6,095 genes). First, we 
determined the average mRNA levels for each of five classes 
of genes, representing DNA deletion, no change, and low-, 
medium-, and high-level amplification (Fig. 4a). For both the 



breast cancer cell lines and tumors, average mRNA levels 
tracked with DNA copy number across all five classes, in a 
statistically significant fashion (P values for pair-wise Student's 
t tests comparing adjacent classes: cell lines, 4 x MM* 1 X 10 -49 , 
5 X 10" 5 , 1 X 10-*; tumors, 1 X 10" 43 , 1 x 10" 2W f 5 X IQ~<\ 
1 X lb~ 4 ). A linear regression of the average log(DNA copy 
number), for each class, against average log(mRNA level) 
demonstrated that on average, a 2-fold change in DNA copy 
number was accompanied by 1,4- and 15-fqldchangesinmRNA 
level for the breast cancer cell lines and tumors, respectively (Fig. 
4a, regression line not shown). Second, we characterized the 
distribution of the 6,095 correlations between DNA copy num- 
ber and mRNA level, each across the 37 tumor samples (Fig. 4b) . 
The distribution of correlations forms a normal-shaped curve, 
but with the peak markedly shifted in the positive direction from 
zero. This shift is statistically significant, as evidenced in a plot 
of observed vs. expected correlations (Fig. 4c), and reflects a 
pervasive global influence of DNA copy number alterations on 
gene expression. Notably, the highest correlations between DNA 
copy number and mRNA level (the right tail of the distribution 
in Fig. 4b) comprise both amplified and deleted genes (data not 
shown). Third, we used a linear regression model to estimate the 
fraction of all variation measured in mRNA levels among the 37 
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tumors that could be attributed to underlying variation in DMA 
copy number. From this analysis, we estimate that, overall, about 
7% of all of the observed variation in mRNA levels can be 
explained directly by variation in copy number of the altered 
genes (Fig. Ad). We can reduce the effects of experimental 
measurement error on this estimate by using only that fraction 
of the data most reliably measured (fluorescence Intensity/ 
background >3); using that data, our estimate of the percent 
variation in mRNA levels directly attributed to variation in gene 
copy number increases to 12% (Fig. 4$Q. This still undoubtedly 
represents a significant underestimate, as the observed variation 
in global gene expression is affected not only by true variation in 
the expression programs of the tumor cells themselves, but also 
by the variable presence of non-tumor cell types within clinical 
samples. 

Discussion 

This genome-wide, array CGH analysis of DMA copy number 
alteration in a series of human breast tumors demonstrates the 
usefulness of defining amphcon boundaries at high resolution 
(gene-by-gene), and quantitatively measuring amplicon shape, to 
assist in locating and identifying candidate oncogenes. By ana- 
lyzing mRNA levels in parallel, we have also discovered that 
changes in DNA copy number have a large, pervasive, direct 
effect on global gene expression patterns in both breast cancer 



cell lines and tumors. Although the DNA microarrays used in our 
analysis may display a bias toward characterized and/or highly 
expressed genes, because we are examining such a large fraction 
of the genome (approximately 20% of all human genes), and 
because, as detailed above, we ate likely underestimating the 
contribution of DNA copy number changes to altered gene 
expression, we believe our findings are likely to be generalizable 
(but would nevertheless still be remarkable if only applicable to 
this set of -6,100 genes). 

In budding yeast, aneuploidy has been shown to result in 
chromosome-wide gene expression biases (13), Two recent 
studies have begun to examine the global relationship between 
DNA copy number and gene expression in cancer cells. In 
agreement with our findings, Phillips et ai (14) have shown that 
with the acquisition of tumorigenictty in an fanraortalized pros- 
tate epithelial cell line, new chromosomal gains and losses 
resulted in a statistically significant respective increase and 
decrease in the average expression level of involved genes. In 
contrast, Platzer et a/. (15) recently reported that in metastatic 
colon tumors only -4% of genes within amplified regions were 
round more highly (>2-fbId) expressed, when compared with 
normal colonic epithelium. This report differs substantially from 
our finding that 62% of highly amplified genes in breast cancer 
exhibit at least 2-fold increased expression. These contrasting 
findings may reflect methodological differences between the 
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studies. For example, the study of Platzer et al (15) may have 
systematically under-measured gene expression changes. In this 
regard it Is remarkable that only 14 transcripts of many thousand 
residing within unamplified chromosomal regions were found to 
exhibit at least 4-fold altered expression in metastatic colon 
cancer. Additionally, their reliance on lower-resolution chromo- 
somal CGH may have resulted in poorly delimiting the bound- 
aries of high-complexity amplicons, effectively overcalling re- 
gions with amplification, Alternatively, the contrasting findings 
for amplified genes may represent real biological differences 
between breast and metastatic colon tumors; resolution of this 
issue will require further studies. 

Our finding that widespread DNA copy number alteration has 
a large, pervasive and direct effect on global gene expression 
patterns in breast cancer has several important implications* 
First, this finding supports a high degree of copy number- 
dependent gene expression in tumors. Second, it suggests that 
most genes are not subject to specific autoregulation or dosage 
compensation. Third, this finding cautions that elevated expres- 
sion of an amplified gene cannot alone be considered strong 
independent evidence of a candidate oncogene's role in tumor- 
igenesis. In our study, fully 62% of highly amplified genes 
demonstrated moderately or highly elevated expression. Tliis 
highlights the importance of high-resolution mapping of ampii- 
con boundaries and shape [to identify the "driving** gene(s) 
within amplicons (16)], on a large number of samplesvln addition 
to functional studies. Fourth, this finding suggests that analyzing 
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the genomic distribution of expressed genes, even within existing 
microarray gene expression data sets, may permit the inference 
of DNA copy number aberration, particularly aneuploidy (where 
gene expression can be averaged across large chromosomal 
regions; see Fig. 3 and supporting information). Fifth, this 
finding implies that a substantial portion of the phenotypic 
uniqueness (and by extension, the heterogeneity in clinical 
behavior) aniong patients' tumors may be traceable to underly- 
ing variation in DNA copy number. Sixth, this finding supports 
a possible role for widespread DNA copy number alteration in 
tumorigenesis .(17, 18), beyond the amplification of specific 
oncogenes, and deletion of specific tumor suppressor genes. 
Widespread DNA copy number alteration, and the concomitant 
widespread imbalance in gene expression, might disrupt critical 
stoichiometric relationships in cell metabolism and physiology 
(e.g., prpteosome, mitotic spindle), possibly promoting further 
chromosomal instability and directly contributing to tumor 
development or progression. Finally, our findings suggest the 
possibility of cancer therapies that exploit specific or global 
imbalances in gene expression in cancer. 
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HER,2/neu Breast Cancer Predictive Testing 

Each year, over 182.0Q0 women in the United States are recurrence risk m women younger than 40 years of a&e for 

diagnosed with breast cancer/and approximately 45.000 .die low- versus Wgh-amplified tumors (54.5% compared to 

of the diseased Incidence appears to be increasing in the 85;7%); this is compared to a recurrence rate of 1 6 7% for 

United States at a rate of roughly 2% per year. The reasons patients whh no HER-2/n eu gene amplification/ HER-2/heu 

for ffie increase ereundear,but non-geneticrisk factors appear stanjif may be particularly importantto establish in women with 

tc^aj^iole^; ^ ^ X? 

-Rve-jfear aundval rates lange^^ppro^tely 65%- ~' Thftttoioe of methofelogy fti determination of HER^2/ 

85%, depending on demographic group, with a significant neu flatus depmis in part on thedinical setting, FDA approval 



I 



percentage, of women experiencing recurrence of their cancer 
within 10 years of diagnosis. One of the factors most predic- 
tive for recurrence once a diagnosis of breast cancer has been 
made is the number of axillary lymph nodes to which tumor 
has metastasized. Most noderpositfre women are given adju- 
vant therapy, which increases their Survival. However, 20%- 
30%*f patients without axillary node involvement also 
develop recurrent disease,, and the difficulty lies in how to i den- 
tify this high-risk subset of patients, These patients could 
benefit -from increased surveillance, early intervention, arid 
treatment * . 

Prognostic markers currently used in breast cancer recur- 
rence prediction include tumor size, histological grade, steroid 
hormone receptor status, DKAploidy, proliferative index, and 
cathepsin D status. Expression of growth factor receptors and 
over-expression of the HER-2/ueu oncogene have also been 
identified as having value regarding treatment regimen and 



1^&$Jx$x (i%> knbwffas c-erbB2) if an oncogene fiat 
encodes a transmembrane glycoprotein that is homologous 
to, but distinct from, the epidermal growth factor receptor. 
Numerous studies. have indicated that high levels of expres- 
sion of this protein are associated with "rapid tumor growth, 
certain forms of therapy resistance, and shorter disease-free 
survival. The gene has been shown to be amplified and/or 
overcxpitssed in 10%- 30% of invasive breast cancers and in 
40%-60% of intraductal breast carcinoma. 3 

There are two distinct FDA-approved methods by. which 
HER-2/neu status can be evaluated: imniunohistochemistry 
GHC, HercepTest™) and FISH (fluorescent in situ hybridiza- 
tion, PathVysion™ Kit). Both methods can be perfoiraed on 
archived and current specimens, the first method allows visual 
assessment of the amount of HER-2/neu protein present on 
die cell membrane. Tbe latter method, allows direct quantifi- 
cation of the level of gene amplification present in the tumor, 
enabling differentiation between low* versus high-amplifica- 
tidh. At least one study has demonstrated a difference in 



for toe Vysis FISH test was granted based on clinical trials 
involving 1549 node-positiye patients. Patients received one 
of three OifFerent treatments consisting of different doses of 
cyclophosphamide, Adriamycin, and 5-fluorouracil (CAF). 
The study showed that patients with amplified HE&»2/neii 
benefited from treatment with higher doses of,adriamycin- 
based therapy, while^those with wimal HBR-2/nea levels did 
not The study, tttersfofe identified a sub-set pf women, who 
because ihey did not benefit from more aggressive treatment, 
did not ne<?d to be exposed to the associated side effects. In 
addition, other evidence indicates that HER-2/neu amplifica- 
tion in node-negative patients can be used as an independent 
prognostic indicator for eafly recurrence, recurrent disease at 
any time and disease-related death* 5 Demonstration of HER* 
2/neu gene, amplification by FISH has also been shown to be 
of value in predicting response to chemotherapy in stage-2 
breast canber patients. - 
. Selection of patieatsf^^^ 

clonal antfl*^1nerapy,^ is based upon demeagca* 
don 6fHEFL-2/neu protein overexpression using HercepTest™ 
Studies using Herceptin c in patients with metastatic breast 
cancer show an increase in time to disease progression, 
increased response rate to chemotherapeutic agents and a small 
increase in overall Survival rate. The FISH assays have not yet 
been approved for this purpose, and studies looking' at response 
to Herceptm e in patients with or without gene amplification 
status determined by FISH are in progress. 

In general, FISH and IHC results correlate well However, 
subsets of tumors are found which show discordant results; 
i.e. t protein overexpression without gene amplification or lack 
of protein overexpression with gene amplification. The Clini- 
cal significance of such resul ts is unclear, Based on the above 
considerations, HER-2/neti testing at SHMC/PAML will utiv 
lize immunohistochemistry (HcrcepTfest 0 ) as a screen, fol- 
lowed by FISH in IHC-negativc cases. Alternatively, either 
method may be ordered individually depending on the clini- 
cal setting or clinician preference. 



AUGUST 1999 



OT code information References 
HEK-2/neuvlaIHC i, 
88342 (including interpretive report) . * 

HEB-2/nen via FISH 3 
88271*2 Mole^larcyto^ 

38274 Molecular cytogenetics, interphase in situ hybrid- 

ization, analyze 25-99 cells 4 

88291. Cytogenetics and molecular cytogenetics, interpre- 
tation and report 

Procedural Information 

Iinmonohistochemistty is performed using the FDA-approVed 
fcAJCG antibody kit^ Hera^tesrG. The DAKO kit contains 
reagents required tp complete a two-step immunohisto* 
cbemi^al staining jmxedqre forjoutinery processed, paraflBn- 
embedded q^imem FoU^^ 
t^rtarrtfbo^tobin^ 

a ready-to-^se dextran-based visoalization reagent This re- 
agent consists of both secondary goat anti-rabbit antibody 
molecules with horseradish peroxidase molecules linked to a 
common dextran polymer backbone, thus eliminating the heed 
. for sequential application of link antibody and peroxidase 
conjugated antibody. Enzymatic conversion of the subse- 
quently added ^hr^mogen results in formation of visible 
reaction product at the antigen site* The specimen is then cotm- 
terstained; a pathologist using light-microscopy interprets 
.results. 

FISH analysis at SHMC/PAML is performed using the 
FDA-^pproved Pa AVywon™ HER-2/neu DNA probe kit, pro- 
ceed by Vysis,lnc> Formalin fixed, paraffin-embedded breast 
tissue is processed using routine histological methods, and then 
slides are treated to allow hybridization of DNA probes to the 
ancle! piesent in the tissue section. thePathvysion^ kit con- 
tains two direct-labeled DMA probes, one specific for the 
tlphridrepetitte 
- ihe chwiir^srjr^^ thefsecondibr the HER- 

^neu oTK^ene located at 17ql U-l 2 (spect^ 
meration of the probes allows a ratio of the number of copies 
of chromosome 17 to the number of copies ofHER-2/ncu to 
be obtained; this enables quantification of low versus high 
amplification levels, and allows an estimate of die percentage 
of cells with HER-2/neu gene amplification. The clinically 
relevant distinction is whether the gene amplification is due 
to incr6«ed gene copy number on the two chromosome 17 
homolbgues normally present or an increase in the number of 
chromosome I7s in the cells. In the majority of cases, ratio 
equivalents less than 2.0 are indicative of a normal/negative 
result, ratios of 2.1 and over indicate that amplification is 
present and to what degree* Interpretation of this data will be 
performed tod reported from theVysis-certified Cytogenet- 
ics laboratory at SHMC 
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